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Abstract 



A class of exact solutions of the field equations with higher derivative terms is 
presented when the matter field is a pressureless null fluid plus a Maxwellian 
static electric component. It is found that the stable solutions are black holes 
in anti de Sitter background. The issue of the stability of the Cauchy horizon 
is discussed. 



I. INTRODUCTION 

Why should one wish to study the internal structure of a black hole? Aside from sheer 
curiosity there is the intellectual challenge of testing our present physical theories up to their 
extreme boundary of validity. 

In fact the spacetime outside a black hole is indeed a rather uninteresting subject. Thanks 
to the no-hair theorem we know the late time structure of the spacetime after the star has 
radiated away all the asphericities [I]. Instead, the ultimate fate of the star that undergoes 
a gravitational collapse, is still an open issue. 

The plausible, still unproven strong cosmic censorship conjecture states @ that the 
singularity at r = is spacelike. The spacetime near the singularity can probably be 
described by the general mixmaster type solution ||. 

It is however puzzling that even an innnitesimally small amount of angular momentum 
generates the Kerr-Newman singularity which is instead timelike. In this case the situation 
is perhaps more dramatic due to the presence of a null hypersurface, the Cauchy horizon 
(CH), boundary of predictability of the evolution of the fields. Like in the Reissner-Nord- 
strom solution, the existence of this null hypersurface makes the existence of a timelike 
singularity in a region beyond it physically irrelevant. 

As it was pointed out by Penrose, |J the Cauchy horizon is also a surface of infinite blue- 
shift. A free-falling observer crossing the Cauchy horizon will measure in-falling radiation 
to have infinite energy density. When only ingoing radiation is present a weak non-scalar 



singularity forms, which is classified as a whimper singularity . Whimper singularities are 
unstable to perturbations that can transform them into stronger, scalar, singularities. 

This scenario can be examined more closely by assuming spherical symmetry. When an 
outflow crossing the Cauchy horizon is included in the analysis, the (Bondi) mass function 
is found to diverge exponentially || at late advanced times v 

m ~ e KV (1) 

where k is the surface gravity of the inner horizon located at v = +00. The only non- 
vanishing component of the Weyl tensor, the invariant ^2, is proportional to the mass 
function and a scalar curvature singularity forms along the Cauchy horizon This 
phenomenon, named "mass inflation", has become a topic of increasing interest and lot of 
effort is devoted to study this subject [§]. 

The Reissner- Nordstrom black hole shares the same causal structure as the Kerr-Newman 
black hole, so it is not surprising that the basic picture derived in spherical symmetry should 
be qualitatively the same for the more general black hole. Some investigations suggest that 
the general scenario derived for spherical symmetry does not change dramatically in a non- 
spherical black hole P-pJjj. In fact one key result of the mass inflation picture is that the 
rate of the exponential divergence of the mass function is entirely determined by the surface 
gravity which is constant along the generators of the Cauchy horizon. 

The relevant question is then to understand how much the classical description of the 
spacetime near the singularity changes when quantum effects are taken into account. In 
particular, would quantum effects enforce the strong censorship conjecture? A first analysis 
as been already performed in |T2[ where the semi classical approach has been used to show 



that the classical picture remains valid up to few Plank lengths from the CH. The mass 
function at that "time" has already grown up to the mass of the observable Universe! 

In fact, the classical description of the black hole interior is simplified by causality: 
behind the event horizon the coordinate r is timelike, so a descent into a black hole is a 
progression in time. The evolution down to any particular radius is only influenced by 
the initial data at larger radii. For this reason, a mean field, semi-classical description of 
quantum effects cannot drastically change the classical picture. This scenario can instead 
change in two-dimensional models, as shown by Jl3 . 



However, there is no reason why one should assume that a perturbative calculation can be 
applied when the curvature has reached planckian level. Only in a framework of a quantum 
gravity calculation one could hope to address this issue. 

In recent years it has become clear that a very convenient way to quantise gravity is to 
describe it as a quantum effective field theory [pL4|— pT6 1 . As a result one is lead to consider 



a more general action that is a functional of any geometrical invariant 3? which can be 
constructed from the principle of general covariance, 3? = {R, R a pR afi C^^C 01 ^ ' , ...}, and 
of the matter field and its derivatives 

Sfag] = J d A x^C{% <p, V 2 K, V 2 0, V 2 ^, V 2 ^,...,) (2) 



2 



This new theory is analogous to the Euler-Heisenberg Lagrangian, low-energy effective the- 
ory of the QED. The important point is that both QED and the Euler-Heisenberg theory 
represent the same theory at different scale lengths. One would not include those new oper- 
ators in the QED Lagrangian since they are suppressed as inverse power of the ultraviolet 
cutoff in the infrared region showing their "irrelevant" character in this case. They can be- 
come instead significant if we move the original cut-off of the theory in the infrared region, 
where the theory must describe phenomena at much lower energy scales. 

This situation is similar in gravity. Due to the smallness of the ratio between the Newton 
constant and the Fermi constant, modifications of Einstein's theory as in eq. (|2|), are experi- 
mentally indistinguishable from the standard theory at ordinary energy scales ||17|| . However 
at higher energy scales it can alter the physics content of the theory in a significant way. 
They might eventually destroy the unitarity of the S-matrix, signalling that a new physics 
sets in at some energy scale. The theory is perfectly well defined below that threshold, and 
it makes precise predictions. This happens also in the so-called renormalizable theories. The 
Standard Model, for instance, has to be regarded as an effective field theory which breaks 
down at the TeV scale, where the Higgs would enter in a strongly self-coupled phase. 

Higher-derivative gravity theories arise also through the coupling of a quantised field 
with the classical background geometry through the renormalization process. They are 
general functions of 3? and are needed to cancel the ultraviolet divergences for any given 
order in perturbation theory. 

However, the standard perturbative approach based on the scaling property of the 
Green's functions under a rescaling of the metric JOJ handles only a finite number of op- 
erators, those which are important for the ultraviolet fixed point. In this way one cannot 
follow the evolution of the irrelevant, i.e. non-renormalizable operators. A coupling is ir- 
relevant, marginal or relevant if it, respectively, it gets smaller, it does not run, it grows 



when the cut-off is lowered from the ultraviolet towards the infrared 20 . The irrelevant 



operators, even if they are not present in the bare Lagrangian, mix their evolution with the 
renormalised trajectory of the relevant couplings. Although near the Gaussian ultraviolet 
fixed point their running is suppressed, they can behave in a quite different manner in other 
scaling region and they can eventually drive the system in a different continuum limit. 

It is therefore important to trace the evolution of all the coupling constants generated 
by the renormalization procedure in order to study the phase structure of the theory. In 
this way one recovers the predictability power of the theory by retaining only few relevant 
operators at a given fixed point. If for example a new scale other than the cut-off is present 
in the problem, the scaling laws change near that scale and the standard ultraviolet relevant 
interactions might not be sufficient to describe the physics at a crossover between ultraviolet 
and infrared. 

The renormalization group approach used in Statistical Mechanics is the best tool to 



study problems where many scales are coupled together |[21||. A powerful way of obtaining 



a differential form of the RG transformation has been formulated by Wegner and Houghton 
Starting from a bare action Sk at the cut-off k one first calculates Sk-Ak in 
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e -s k -AkW = J p^] e -^[WI ( 3 ) 

by using the loop expansion, where ip and respectively have non-zero Fourier components 
only in the momentum shells k — Ak < p < k and p < k — Ak. The differential RG 
transformation is obtained by taking the limit of an infinitesimal shell Ak/k —>■ bkjk. The 
higher loop contributions in Eq. (|3]) are suppressed as powers of Ak/k in the limit Ak — > 
for finite k and an exact, non-perturbative, one-loop RG equation is obtained 

dk ~ 2 { 5^ ) + { 5<P y 5^ ) 5<p > [ ] 

where the brackets indicates the sum over the Fourier components within the shell. This 
functional equation rules the evolution of all the interaction terms that are generated in the 



renormalization procedure. It has been applied in [23| to discuss the evolution of higher- 
derivative (HD) operators generated by the inflation field. In particular it has been shown 
that even if they are not present at the cut-off - which has been chosen below the Planck 
scale! - they show up above the mass threshold and influence the renormalised flow. 

A consistent quantum gravity calculation in the framework of the effective theories, must 



be performed by including the HD operators from the beginning |L6[ . Within a perturbative 
(weak field) calculation one can use the Einstein-Hilbert truncation [24], but in a strongly 
non perturbative situation one should consider a more general Lagrangian and, with the 
help of the renormalization group equations, look for the consistent field configurations that 
effectively dominate the path integral at that scale of energy. 

The first step towards this goal is to understand the tree-level structure of the space- 
time in the interior of a black hole, when higher-derivative operators are included from the 
beginning. 

Those new operators can produce new istantons configuration in the path integral formu- 
lation of quantum gravity, with stronger weight than the solutions of the standard Einstein 
theory. If one wants to perform a saddle point evaluation of the path integral the first step 
is to determine the stationary points of the (Euclidean) action. Expansion around those 
saddle points are performed by writing 

liw = 5V + h/iv (5) 

where, as usual, treating g^ u as a background field and h^ v as a quantum field one generates 
the usual perturbation expansion that can be expressed in terms of Feynman diagrams. The 
strong assumption that is usually believed, is that vacuum state for the renormalised system 
is still given by the g^ v field configuration. Is this assumption still valid in the renormalised 
system when new higher dimensional operators are generated by the renormalization pro- 
cess? If not, we have chosen a wrong, unstable, saddle point. Strominger has shown |25j 
that for quadratic gravity, the Minkowsky solution is the only one having zero ADM energy. 
However that is not necessarily true for a more general non-linear gravity theory. As it has 
recently been pointed out by p6| where a a completely new vacuum structure has appeared 
already for R + R 3 a Lagrangian. 
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In this work we shall reduce the original HD theory to the standard second order form by 
means of the method outlined in p6|j27[| . We shall obtain an equivalent (locally isomorphic) 



second order theory with an additional, non-geometric degree of freedom represented by a 
non-ghostlike scalar field. We call "vacuum" a solution which is stable with respect to the 
elementary excitation of this new field. We show that for a general non-linear Lagrangian 
the theory has non-trivial vacuum solutions representing BH embedded in Anti de Sitter 
spacetime. 

We shall also comment on the stability of the Cauchy horizon in those solutions. 



II. BASIC EQUATIONS 

Let us consider the following action 

S = J d 4 x^g-(l g (R) + l m (A, g)) (6) 

where l g (R) is a non-linear function of the scalar curvature R which we shall suppose to be 
an analytic function of R 

l 9 (R) = ±^R n (7) 

e n are the coupling constants of the generic power of R and l m {A) is the Lagrangian density 
of the matter field A. We use the signature (-+++) and we set c = Ti = 8nG. With these 
units e = —A corresponds to the vacuum energy, being A the cosmological constant and 
ei = 1/2. The gravitational field equations for such a system are the following fourth order 
differential equations 

^g~l g (R)} = l{R)R pv - h g (R)g, u - V.VJ'JR) + g,MR) = T^A, g) (8) 



where is the stress-energy tensor of the matter field which obeys V = 0. 



In order to analyse the above theory it is convenient |26| to introduce an auxiliary field 
ip so that the action now reads 

S = Jd 4 xV=g(i' g (R)(R-ii)-iM+UAg)) (9) 

with the equation of motion 

JJ(V)OR-V) = o. (io) 

The new action, written with the help of the auxiliary field, is (at the classical level) equiv- 
alent to the original theory provided that we are not at the critical points l g (R) = 0, but 
it is not of the second order form. In order to reduce the action ( |10|) in a canonical second 
order form, one introduces a pair of new variables (g^ u ,uj) related to g^ u by 
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(11) 



where I' > is assumed and ^f(u) becomes a solution (not necessarily unique), of 



WW) - e U 



0. 



(12) 



This transformation has already been used in the literature, and it has been generalised to 
the case of Lagrangians depending on \7 2k R (see [27 for a general review on the subject). 
Under the above conformal transformation the Ricci scalar density transforms like 



^g-R = y/=7je-"(R - |(Vu;) 2 - 3V 



'CO 



(13) 



By dropping a total derivative term, the action has been reduced to the canonical form 

S = J d 4 Xv ^(R - ^ v d^d v LO - 2V(lu) + e- 2uJ l m (A, e^g)) (14) 
where the potential term 

V(u) = ^(e^(uj)-l g ^(uj))). (15) 
depends on the specific original higher-derivative theory. The field equations thus read 



G 



t + T 



where Gnu is the Einstein tensor of the reduced theory 



G Rfii/ 



tn U is the stress energy-tensor of the scalar field 



3r 



d^ujd v uj - -g^g^daudpuj - g^ v V{uj) 



(16) 



(17) 



(18) 



and T^y an effective stress-energy tensor generated by the interaction between the original 
matter field and the non-geometric gravitational new degrees of freedom represented by the 
scalar lu, 



(19) 



One should also note that in general the two terms on the left hand side of the equation 
(|i~6|) are not separately conserved, but it holds instead the conservation law V JG^ = 0. We 
stress that the spin-0 field introduced with the help of eq.([TT|) has physical meaning since 
it corresponds to an additional degree of freedom already present in the starting non-linear 
Lagrangian. Therefore the advantage of having used this approach is that this new degree 
of freedom has been now made explicit in the reduced Lagrangian. 
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Since we are interested in studying the structure of the CH in the presence of HD terms, 
we consider a system where the gravitational field is coupled with a pressureless null fluid 
plus an electromagnetic contribution coming from a static electric field generated by a charge 
of strength e. The stress-tensor for the matter field is therefore given by 

Tp, = plpl v + (20) 

where E v = e 2 /87rr 4 diag(l, 1, —1, —1) is the Maxwellian component of the electric field, p is 
the energy density of the radiation and I ^ is a null vector tangent to the radial null geodesic. 
In the following we shall specify our calculations in the case of spherical symmetry. One can 
introduce a null tetrad n M , m^m^} with 1^ = —1 = —m^rh^. The metric tensor then 
reads 



9vv 

and 



m^m u ) - l( p n v) . (21) 



e 2 



E »v = -]^lhv n ») + m ^ m v)\ (22) 

In our case the matter Lagrangian is conformally invariant, and the equation of motion 
for the matter field completely decouple from those for the scalar field. We analyse the 
case of us = const, and the solutions of the equation of motion correspond to constant field 
configuration that are extrema for the potential 

dV , , 



The vacuum case -no matter field- has been analysed in ||26|| . In our case we see that the 
conformal transformation in ([LTD generates a conformal rescaling of the null tetrad in (j21~D 
with 1^ — > Z M = exp(u;/2)Z At . It is then easily seen that the stress-energy tensor in the 
reduced theory in eq.(|l9D has the same functional form of that in the original theory, with 
the "renormalised" energy density and charge 

p — > p = exp(—2u)p e 2 — ► e 2 = exp(— 2u;)e 2 . (24) 

The solutions of the new field equations (|T|) can therefore be classified by looking at the 
minima of the potential term (15). Even if in the original theory a cosmological term is 



not present, the solutions of the reduced theory are given by a Vaydia like spacetime in a 
de Sitter (dS) or Anti-de Sitter (AdS) background, depending on the value of the potential 
at minimum. In particular, one can use the Eddington-Finkelstein coordinates so that an 
explicit form of the metric element in the reduced theory reads 

ds 2 = dv{2dr - fdv) + r 2 dtt 2 . (25) 

where 
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2m(v) e 2 A o 

f = 1 -—tr 1 + -2-3 r (26) 



is the so called mass function, A = V(u c ) and u c is the location of any extrema of the 
potential. In a local chart adapted to the inner horizon the Cauchy horizon is located at 
v = +00, and / = 0. If we compare this metric with the solution we would have gotten in 
the Einstein gravity, we note that, according to eq. (|24]) the charge and the energy densities 
are screened or anti-screened depending on the sign of u c . It should be also noted that, 
while in standard gravity one has always AdS solution of the kind (j25D in our case the above 
metric is a solution of the field equations only for special values of the cosmological constant. 

The global spacetime structure of those solution is well known. An updated review can 
be found in the contribution by Chris Chambers in these proceedings, while recent results 



on topological (AdS) black holes can be found in the lecture of R. Mann |38 |. See also there 
for the figures displaying the conformal structure of the two spacetimes. 

Here we shall simply stress that the black hole spacetime in both dS and AdS solutions 
may have Cauchy horizons depending on the roots of the equation 

/(r) = (27) 

If the dS case is considered eq.(p7|) may have three positive roots and the spacetime has 
therefore a cosmological horizon besides an inner and an event horizon. Thanks to the work 



of Brady and Poisson [39] we know the Cauchy horizon is stable if the surface gravity of the 



cosmological horizon is greater than the surface gravity of the inner horizon 

K C S > k C h (28) 

The surface gravity of the horizon in the original theory is obtained by a simple scale trans- 
formation (from the definition of surface gravity), and one can conclude that the stability 
criteria (PS"D holds in the original theory as well. 

However, a black hole in dS background does not represent a vacuum solution since the 
field sets on a maximum of the potential. The stable vacuum solutions for the excitations 
of the u field are instead given by black holes in AdS background and for those solutions 
the Cauchy horizon is a subtle question. The spatial infinity is now timelike, the radiative 
falloff at late times is not well understood []30|| . For a power law decay of the kind 



5( v ) = l/(av) p (29) 

with p > 3 and a has to be introduced on dimensional grounds, the Cauchy horizon is 
unstable. In fact in general the energy density p b s measured by an observer that crosses 
the Cauchy horizon would diverge as 

p obs ~ 8 exp(2K C Hv) (30) 

thus signalling the instability of the Cauchy horizon. The CH in the original theory would 
also be unstable since the conformal transformation in eq. (|ll]) is regular. However for an 
exponential decay rate of the kind 
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5 ~ exp(— 2atv) 



(31) 



one would not see any divergence for 

a > k C h (32) 

In a more realistic situation, one should use the eq. (|29| ) to mimic the radiation falling 
into the event horizon. Outside the event horizon there is no reason to suppose that the 
background being AdS. Our effective "cosmological" constant is therefore dynamically gen- 
erated during the gravitational collapse for the presence of the quantum fluctuations when 
the curvature has reached planckian level. The initial data coming from larger radii is still 
given by the classical evolution and therefore the Cauchy horizon would still be unstable 
even in the presence of a negative (or positive) A term in the metric. 



III. R 2 , R 3 + R 2 GRAVITY 

Let us analyse some cases in detail. If we consider R 2 gravity, the density Lagrangian of 
the gravitational field reads 

l g (R) = e + e x R + ^e 2 R 2 . (33) 

then, eq. (p"2|) can be inverted to yield ip = ln(l/2 + e 2 '0) provided ip > — l/2e 2 . In particular 
the potential term becomes 

^) = ^(l-e ie -) 2 -e- 2 % (34) 

We see that for a positive cosmological constant in the original theory the potential term is 
not bounded from below and there are no stable solutions of the quadratic gravity theory. In 
the case of a negative cosmological constant, the potential term has only one minimum for 
uj c = — In 2 and we find the BH solutions given by eq. fl2"6"D . We note that the " renormalised" 
charge and energy density have now an increased strength of a factor 4 according to (p4[). 

Another interesting case is that of R 3 gravity where we suppose that the original La- 
grangian has the following functional form 

l g (R) = e 1 R+^e 2 R 2 + ^e 3 R 3 (35) 



in this case we can invert the relation (12) in order to obtain the two solutions 



^±(w) = -( - e 2 ± y/e 2 2 + 2(e» - ei )e 3 )) (36) 

in the two branches ip > — e 2 /e^ and ip < — e 2 /e 3 and for the reality condition we also 
consider uj > oj* = ln(l/2 — e^/es). The potential term is therefore given by 
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V{^)) = |^Me- 2 ^(l + g|^(o;)) (37) 

In the first case the potential V(if)+(u) has one minimum for u c = — In 2 with V(o; c ) = 
and a local maximum at another uj c < with V(u; c ) > 0. Therefore in the ip + branch we 
can conclude that the solution of the field equations is given by a Vaidya like metric in a 
dS background. In the ip- branch we have only one local minimum that corresponds to a 
negative value of the potential, and therefore the metric is still a Vaidya like, but in AdS 
background. 

Other interesting cases are given by Lagrangian like l g (R) = R + bR\n(a + R), with a 
and b constants, which would arise from a one-loop contribution in an interacting stress- 
energy tensor. In any case one can use the above outlined procedure, and derive the explicit 
expression of the dS or AdS black hole solution. 

The use of the conformal transformation is not the only method to obtain the solution 
of the original theory. There is also a direct procedure which is of course equivalent. Let us 
consider the field equations given in @ and set 

R„ v = C(R) ( P y v + E^) + -Rg^ (38) 

with constant Ricci scalar R. The following statement is true: 

The solutions of the field equations in ([<|) with the stress energy-tensor in 
Ricci tensor of the form $3$) have 

Rl' g (R) - 2l g (R) = (39) 
l' g (R)C(R) = 1. (40) 

The above conditions have to be considered as implicit relations that constrain the possible 
values of the constant R. They can be deduced by direct substitution of fl38D in @ and 
by projection along l^n u and n^n u respectively. One also notes that the first of the above 
relations in ( |38|) is equivalent, aside from a constant factor, to the extrema condition for the 
potential, and the second relation can be read off as C(R) = e~ u which is the first of the 
relations in (|ll|). This shows the equivalence of the two approaches in obtaining solutions 
of the HD theory. However the conformal transformation method is simpler and elegant, 
while the direct method has the advantage that can also be used when lg(R) = 0. The set of 
consistency equations in (|38D can also be deduced by considering the following stress energy 
tensor 

T^ = C(R)- 1 pl li l v + E ia , (41) 
and with the Ricci tensor given by 

R„u = pl^lu + C(R)E llu + -Rg^. (42) 



2C ) and the 
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The above relations can be used to investigate the full global content of the theory. Once 
eq.fl38|) has been solved, one introduces locally the conformal frame where the new degrees 
of freedom is made explicit. Then one can address the question of local stability and mass 
of the field. 



IV. CONCLUSIONS 



As we have stressed in the introduction, the non-linear modifications of the standard 
Einstein theory are needed for taking into account of the quantum effects when the Weyl 
curvature has reached Planckian levels. We have seen that stable black hole solutions appears 
in these cases. What is the possible role of these new solutions in the standard mass inflation 
scenario? It is interesting to understand the structure of the CH singularity in a more 
consistent calculation where these new terms are dynamically generated during the collapse. 
We have argued that the role of those new terms would not change the (unstable) character 
of the Cauchy horizon. We therefore think that the mass inflation phenomenon also takes 
place in such a background. However since the appearance of an effective A term changes 
both the location of the horizon and the surface gravity of the inner horizon, a detailed 
calculation must be performed in order to address this issue. 
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